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There was interest in understanding the impact of out-of-round nozzle extension on the 
nozzle side load during transient startup operations. The out-of-round nozzle extension 
could be the result of asymmetric internal stresses, deformation induced by previous tests, 
and asymmetric loads induced by hardware attached to the nozzle. The objective of this 
study was therefore to computationally investigate the effect of out-of-round nozzle 
extension on the nozzle side loads during an engine startup transient. The rocket engine 
studied encompasses a regeneratively cooled chamber and nozzle, along with a film cooled 
nozzle extension. The computational methodology is based on an unstructured -grid, 
pressure-based computational fluid dynamics formulation, and transient inlet boundary 
flow properties derived from an engine system simulation. Six three-dimensional cases were 
performed with the out-of-roundness achieved by three different degrees of ovalization, 
elongated on lateral y and z axes: one slightly out-of-round, one more out-of-round, and one 
significantly out-of-round. The results show that the separation line jump was the primary 
source of the peak side loads. Comparing to the peak side load of the perfectly round nozzle, 
the peak side loads increased for the slightly and more ovalized nozzle extensions, and either 
increased or decreased for the two significantly ovalized nozzle extensions. A theory based 
on the counteraction of the flow destabilizing effect of an exacerbated asymmetrical flow 
caused by a lower degree of ovalization, and the flow stabilizing effect of a more symmetrical 
flow, created also by ovalization, is presented to explain the observations obtained in this 
effort. 
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Nomenclature 


Ci,C 2 ,Cs,C f i= turbulence modeling constants, 1.15, 1.9, 0.25, and 0.09. 

C p = heat capacity 

D = diffusivity 

F yZi F y> F z = integrated force, and component forces in the lateral direction 

/ = frequency 

FI = total enthalpy 

K = thermal conductivity 

k = turbulent kinetic energy 

L/S = ratio of long-axis to short-axis 

Q = heat flux 

T = temperature 

t = time, s 

u = mean velocities 

V 2 = Z u 2 

x = Cartesian coordinates or nondimensional distance 

a = species mass fraction 

e = turbulent kinetic energy dissipation rate 

6 = energy dissipation contribution 

= viscosity 

= turbulent eddy viscosity (=pC M k 2 /s) 

FI = turbulent kinetic energy production 

p = density 

a = turbulence modeling constants, 0.9, 0.9, 0.89, and 1.15 for Eqs. (2), (4), (5), (6). 

r = shear stress 

co = chemical species production rate 

Subscripts 
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r 


= radiation 


t = turbulent flow 


I. Introduction 

Nozzle lateral forces during startup and shutdown transients, if not properly managed, are known to cause 
severe structural damages to the engine hardware to almost all liquid rocket engines during their initial development 
[1-4]. Therefore, transient nozzle side load is considered a high risk item and a critical design issue during any new 
engine development. For that reason, many research efforts [5-18] have been devoted to understanding the side load 
physics and their impact on the magnitude of side loads. For regeneratively-cooled engines such as the Space Shuttle 
Main Engine (SSME), the peak side load generating physics have been identified as the X shock oscillation across 
the nozzle lip [7]. For film-cooled engines such as the Japanese LE-7A engine and the U.S. J-2X engine, the major 
side load generating physics have been associated with the jump of the separation line [3, 8]. Other side load physics 
such as the Free-Shock Separation (FSS)-to-Restricted-Shock Separation (RSS) transition have been mentioned as 
the critical one for the European Vulcain engine [9]. 

With respect to the J-2X engine under development, the impact of the side loads caused by the jump of the 
separation line has been considered. However, liquid rocket engine nozzles, being large, relatively light weight 
structures, are probably never truly round. The effect of nozzle out-of-roundness, on the transient nozzle side loads 
is therefore a potential issue and needs to be explored. The causes of nozzle out-of-roundness could be, but are not 
limited to, asymmetric material internal stresses induced in previous tests, the asymmetric loads induced by other 
hardware attached to the nozzle, and nozzle wall material deformation, such as creep, induced in previous engine 
tests. 

Our first attempt on this subject, report in Ref. [19], studied the effect of out-of-roundness over the entire thrust 
chamber (main combustion chamber, nozzle, and nozzle extension), with four different degrees of ovalization: a 
perfectly round, a slightly ovalized, a more ovalized, and a significantly ovalized. All three ovalized cases were 
elongated along one of the lateral coordinate, the z-axis only, since it was postulated that the elongation along the 
other lateral coordinate y-axis is unlikely, due to the attached structure of the turbine exhaust manifold (TEM). The 
results show that the separation line jump was source of the peak side load for the first three configurations, with the 
peak side load increasing as the degree of out-of-roundness increased. For the significantly out-of-round case, the 
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peak side load was reduced to a level comparable to that of the round nozzle due to a mechanism named splitting of 
the separation line jump [19]. 

In this effort, we have shifted our attention from ovalizing the entire thrust chamber to that of the nozzle 
extension only. Nozzle extension, being structurally thin and not containing cooling tubes, is probably more likely to 
deform than the other components such as the regeneratively-cooled main combustion chamber and nozzle. Since 
the elongations in both the lateral y and z axes are now likely, six three-dimensional computations were performed 
with the out-of-roundness achieved by three different degrees of ovalization, and elongated on both y and z axes: 
one slightly out-of-round, one more out-of-round, and one significantly out-of-round. In order to compare results 
with those of previous study, the maximum deformations of the three different degrees of ovalization are identical to 
those used in the previous study [19]. The results of these computations are presented and discussed herein. 


II. Computational Methodology 

A. Computational Fluid Dynamics 

The computational fluid dynamics (CFD) methodology is based on a multi-dimensional, finite-volume, viscous, 
chemically reacting, unstructured grid, and pressure-based formulation. Time-varying transport equations of 
continuity, species continuity , momentum, total enthalpy, turbulent kinetic energy, and turbulent kinetic energy 
dissipation were solved using a time-marching sub-iteration scheme and are written as: 
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A predictor and corrector solution algorithm was employed to provide coupling of the governing equations. A 
second-order central-difference scheme was employed to discretize the diffusion fluxes and source terms. For the 
convective terms, a second-order upwind total variation diminishing difference scheme was used. To enhance the 
temporal accuracy, a second-order backward difference scheme was employed to discretize the temporal terms. 
Point-implicit method was used to solve the chemical species source terms. Sub-iterations within a time step were 
used for driving the system of second-order time-accurate equations to convergence. Details of the numerical 
algorithm can be found in Refs [20-23]. 

An extended k-s turbulence model [24] was used to describe the turbulence. A modified wall function approach 
was employed to provide wall boundary layer solutions that are less sensitive to the near-wall grid spacing. 
Consequently, the model has combined the advantages of both the integrated-to-the-wall approach and the 
conventional law-of-the-wall approach by incorporating a complete velocity profile and a universal temperature 
profile [25]. A 7-species, 9-reaction detailed mechanism [25] was used to describe the finite-rate, hydrogen/oxygen 
afterburning combustion kinetics. The seven species are H 2 , 0 2 , H 2 0, O, H, OH, and N 2 . The thermodynamic 
properties of the individual species are functions of temperature. The multiphysics pertinent to this study have been 
anchored in earlier efforts, e.g., SSME axial force and wall heat transfer [20], SSME startup side load and dominant 
shock breathing frequency [7], J-2X startup and shutdown side loads for a nozzlette configuration [8], nozzle film 
cooling applications [26], and conjugate heat transfer [27], 

B. Transient Startup Sequences 

Transient system-level simulation is a vital part of the computational methodology, because it provides the time- 
histories of the inflow properties entering the nozzle. Simply put, the ramp rates, or histories, of the inlet pressure, 
fluid temperature and species concentrations play an important role in determining the type of side load physics, 
magnitude and duration of the side loads during the transient operations. In terms of the type of side load physics, it 
was demonstrated that a combustion wave occurred during a transient startup of the SSME due to the excess fuel 
dump [7] ; On the other hand, it was also demonstrated that the combustion wave disappeared when the excess fuel 
dump was replaced with helium purge in the J-2X startup transient [8]. In terms of the magnitude and duration of 
the side loads, it was demonstrated that the peak side load dropped to about half after the total ramp time decreased 
from 5 s to 1 s [7], for the SSME startup simulation. In addition, the combustion wave captured in the 5 -second 
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nominal ramp time was eliminated when the ramp time was shortened to 1 s [7]. In other words, the time-varying 
inlet flow properties determine the all-important residence time of the side load generating flow inside the nozzle. 

The system-level simulation is based on a lumped, control-volume analysis approach to model the rocket engine 
as a network of components and sub-components. This method of transient system-level modeling has been shown 
to be effective in simulating the low-frequency, transient physics associated with the operation of previous and 
existing rocket engines (SSME, RL-10, IPD, etc.) and therefore, is an important tool in the design and planning of 
sequencing the transient events of rocket engine operation. Figure 1 show the engine system transient model 
computed inlet fluid pressure and temperature time-histories, for the main combustion chamber (MCC) and the 
turbine exhaust gas (TEG) flows, during a J-2X startup transient. TEG flow is used as film coolant for the J-2X 
nozzle extension. The transient reactant species composition (not shown) obtained from system modeling was 
preprocessed with the Chemical Equilibrium Calculation program [28], assuming the propellants were ignited to 
reach equilibrium composition immediately downstream of the injector faceplate and turbine exhaust manifold inlet 
boundary, for the respective MCC and TEG flows. This J-2X engine start transient process involves thermal-fluid 
physics phenomena and safety-based operating practices that are typical of a conventional LH2/LOX rocket engine 
based on a gas generator (GG) power cycle with spin-start assistance for the turbomachinery. 

During engine start, the GG fuel and oxidizer injectors are both flowed with helium which serves as a purge or 
buffer flow. This buffer flow reduces the possibility of backflowing of either propellant into the other propellant’s 
manifold during the process of establishing the propellant flows to the GG. After initiation of the GG injector 
purges, the main fuel valve is the first primary engine valve opened, establishing the initial fuel flow in the main fuel 
injector and MCC augmented spark igniter (ASI), and thus ensuring a fuel-rich start of the MCC and safely allowing 
the J-2X fuel turbopump to start increasing in power. The helium spin-start system is then activated, providing the 
initial break-away momentum to respective pump shafts and causing the increase of pump flow-rates, by injecting a 
short-duration (~1.4 sec), high-pressure/high-flow of helium gas into the inlet of the fuel turbine, immediately 
downstream of the GG discharge. The spin-start helium flow back-pressures the GG injector manifolds prior to 
propellant flowing to the GG and thus delaying the start of the GG. Simultaneous with the initiation of the helium- 
spin assist system, the main oxidizer valve is opened to allow priming of the downstream oxidizer system with 
liquid oxygen. Also in this same period, the MCC ASI oxidizer valve is opened, initiating low-flow/low-pressure 
combustion gas in the MCC ASI chamber and MCC during this early timeframe. In very rapid manner, as the spin- 
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start flow diminishes after the helium spin-start valve is closed, the back-pressure to the injectors decays, and the 
GG injectors are relieved to flow. By precise sequencing of the GG valves in coordination with the helium spin- 
start valve and the main propellant valves particularly, a fuel-rich lead flow to the GG is established followed shortly 
after by oxidizer flow to the GG. Once oxidizer and fuel flows are established and within a specified mixture ratio 
range, the GG pyrotechnic igniter is sequenced to pulse and start combustion of the GG flows. Thus, the processes 
of the downstream injection of high-pressure helium spin-start flow, GG ignition, the ramp-up of pressure and flow 
rate throughout the engine, and the priming of 
the main oxidizer injector with high-density 
oxygen (supercritical, but typically referenced 
as "liquid” when describing the injector- 
priming process) occur in a nearly 
simultaneous manner but yet in a very 
coordinated, chronological sequence that is 
required for a safe engine start transient. It is 
during the spin-start system operation, prior to 
any significant flow through the MCC main 
injector, that the nozzle side-loads hit a peak 
magnitude (at 0.80-1.10 sec). 

Further downstream, during the activation 
of the GG injectors’ helium purges and prior 
to spin-start flow initiation, a very low 
pressure helium flowfield is established in the nozzle and MCC in the otherwise ambient vacuum environment. This 
initial flowfield is established via the TEG flowpath, which discharges into the nozzle at the junction of the regen 
nozzle and nozzle extension. This helium flowfield is further effected as the helium spin-start flow is activated and 
exhausts into the MCC-nozzle flowfield via the TEG flowpath (after passing through the fuel and oxidizer turbines). 
The helium in the MCC-nozzle flowfield is later pushed out after the decay of the spin-start flow as the initial fuel 
and low-flow/low-pressure igniter flows become established to the MCC. However, the discharge of helium spin- 
start flow from the TEG into the nozzle flow field occurs prior to any significant flow through the main injector. 



Fig. 1 Computed inlet pressure and temperature histories for 
the main combustion chamber and turbine exhaust gas flows 
during the start-up transient. 
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thus promoting the possibility of loading the nozzle asymmetrically. As the spin-start helium mass flow decays, the 
GG fuel and oxidizer valves are opened, causing build-up of the GG hot gas flow rate, thus further accelerating the 
pumps and powering the buildup of both fuel and oxidizer pressure and flows throughout the engine. This transition 
from helium spin-start gas to GG-generated hot gas is not instantaneous, and there is a short period of time during 
which engine system flows and pressures decay. 

During the turbine-flow transition period, the oxidizer dome of the main injector primes with liquid oxygen, 
which is a crucial step in starting the engine safely, as the oxidizer flow entering the MCC then becomes more 
controllable and predictable, increasing the 
likelihood for stable MCC combustion as engine 
power is sharply increased. As the flow through 
the main oxidizer injector is initiated and then 
primes to a fully liquid flow, the MCC pressure 
rises sharply to the first small plateau around the 
1.4- 1.6 sec timeframe as shown in the transient 
MCC pressure profile in Fig. 1. The rise in MCC 
pressure correspondingly causes a decrease in 
overall fuel system flow and thus decreased fuel 
to the MCC as well, thus causing the brief 
stagnation or plateau in MCC pressure. Figure 1 
shows that the majority of the increase in MCC 
pressure and temperature occurs between 1 .4 and 
3 s. 

III. Computational Grid Generation 

Figure 2 shows a grid layout of the J-2X engine computational model. It includes the MCC, the nozzle, the 
turbine exhaust manifold (TEM), and the nozzle extension. The TEM consists of an inlet duct and a torus in which 
the TEG flow splits, flowing around the torus and subsequently into the nozzle extension. The exit of the TEM is a 
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narrow annular slot of constant area through which the TEG exhausts onto the nozzle extension inner wall [19, 30], 
A separate steady-state mainstage analysis with this TEM configuration indicated it produces a 9% higher, than 
uniform, mass flow into the nozzle extension in the region under the TEM inlet duct, for the round nozzle. The 
pressure monitoring ring shown in Fig. 2 indicates the locations of the monitoring points, where the circumferential 
pressure distribution was monitored. 

The three-dimensional thruster grid, shown in Fig. 2, was generated by rotating an axisymmetric grid first 
without the TEM. The TEM grid was constructed separately as it is asymmetric to the central axis. The grid was 
completed by merging these two grids at an interface. The software package GRIDGEN [29] was used to perform 
the grid generation. The general layout of the outer boundaries and the wall boundaries of the MCC, nozzle, TEM, 
and nozzle extension are similar to those of a previous J-2X side load study and are, described in detail in Ref. [8], 
The rotated grid ensures the grid is symmetrical to the central axis and therefore the computed lateral forces came 
from the asymmetrical flow physics and not from asymmetrical grid cells. 

A grid study procedure was developed for transient nozzle side load computations following the results of a 
SSME transient nozzle side load study [7], by assuming the grid study can be performed on steady-state axial force 
computations. This is necessitated by the cost-prohibitive computational requirement for computing a three- 
dimensional, transient, turbulent reacting flow for 3 to 5 s, especially for sea level simulations since both the flow 
residence time and chemical reaction time are long. In addition, since the peak side force is usually a small 
percentage of the axial force, it is further assumed that if the grid density is suitable for the steady-state axial force 
calculations, then it is likely adequate for the transient side force analyses. This procedure was demonstrated 
adequate as described in the SSME side load study [7] where a total of 1,275,120 grid points and a circumferential 
division number of 72 were used. The predicted peak side loads and dominant frequency for sea level startup 
transient agreed very well with those of tests [7], as summarized in Table 1. The two local peak side loads of 80, and 
102 kN were computed from FSS-to-RSS and RSS-to-FSS transitions, respectively. The two dominant frequencies 
of 122 and 125 Hz were reduced from pressure and heat flux signals, respectively. 


Table 1 A comparison of past CFD results with those of engine tests 



Engine 

CFD 

Engine 

Test 

Physics 

Ref. 

Local Peak F yZj 

SSME 

212 

SSME 

200 

X shock oscillation over lip 

7 


9 




kN 

SSME 

80, 102 

SSME 

90 

shock transitions 

7 

Dominant/, Hz 

SSME 

122, 125 

SSME 

120 

X shock oscillation over lip 

7 

Local Peak F yz 

J-2X 

249 

LE-7A 

259 

separation line jump 

8 

kN 

J-2X 

80 

LE-7A 

102 

FSS-to-RSS transition 

8 


In our first J-2X side load study [8], three grids with 1,581,306 points, 2,011,912 points, and 2,058,192 points 
were generated, all with circumferential division number of 72. The three grids were tested for axial force 
predictions and that of 2,058,192 points was selected. The computed local peak side loads for startup transient at sea 
level compared reasonably well with 


those of a similar engine, also 
summarized in Table 1. Note that 
unlike the SSME which is a 
regeneratively-cooled engine, both 
the I-2X and LE-7A engines are 
film-cooled engines. As a result, the 
peak side load physics switched from 
X shock oscillation over the nozzle 
lip for the SSME, to that of 
separation line jump for the J-2X and 
LE-7A engines. Subsequently, due to 
the availability of more 
computational resources, a grid size 
of 4,421,166 points with a 
circumferential division number of 
120 was used in Ref. [19], and a grid 
size of 20,222,021 points with a 
circumferential division number of 



Fig. 3 Cross-sectional view of the six ovalized thrusters. 


240 was used in Ref. [30]. In this study, a grid size of 4,421,166 points with a circumferential division number of 
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120 was selected for all the cases, due to its better computational efficiency than that of the 20,222,021 points grid, 
and its better resolution in teepees capturing [30] than that of the 2,058,192 points grid. 

In this study, the out-of-roundness was achieved by ovalizing the nozzle extension. Figure 3 shows the cross- 
sectional view of the surface of the six ovalized nozzle extensions. Cases 1, 3, and 5 are the slightly ovalized, more 
ovalized, and significantly ovalized nozzle extensions that were elongated in the z-direction, respectively; while 
Cases 2, 4, 6 are those elongated in the y-directions. For convenience of discussion of the results, the x- and y- axes 
are placed on each of the cross-sectional view of the cases, with the origin placed at the geometrical center. The 
TEG inlet duct is placed on the y-axis and its azimuthal angle count is zero degree. The angle goes clockwise from 
zero to 180 degrees, and goes counter-clockwise from zero to -180 degrees for the circumferential pressure profiles. 
It is noted that although the cross-sectional flow areas of the four nozzles do not have to be the same, it is further 
assumed that those are the same such that the results are compared on an equal flow area basis. The quantitative out- 
of-roundness follows those of a previous study [19] in which the entire thrust chamber was ovalized, hence the 
results can be compared on an equal degree of out-of-roundness basis. That is, ratios of the long axis to short axis 
(L/S) are 1.0086, 1.0346, and 1.4400. The slightly out-of-round case, with a maximum deformation of ±0.25 in, is 
the most likely scenario out of the three. The more out-of-round case, with a maximum deformation of ±1.00 in, is 
also a likely scenario. The significantly out-of-round case, with a maximum deformation of ±11.6 in, is a 
hypothetical case. The purpose of the significantly out-of-round case is to exaggerate the severity of the deformation 
in order to capture the hidden physics that otherwise may be hidden in the other two cases. 


IV. Boundary and Inlet Conditions 

Since J-2X is the upper stage engine of the Ares vehicles, fixed freestream boundary conditions were set 
corresponding to 100,000 ft. Time-varying inlet flow boundary conditions were used at the MCC and TEG inlets. 
These inlet flow properties (obtained from the system level simulations) were the time varying total pressure, 
temperature and reactant compositions. For engine startup computations, the thermal wall boundary condition was 
initially adiabatic. When the startup transient simulation reached 1.4 s, the point where the MCC pressure started to 
ramp up, wall temperature profiles obtained from a separate steady-state calculations were imposed onto the MCC 
and nozzle wall. The nozzle extension and TEM walls remain adiabatic (discussed in [7, 8]). 
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V. Results and Discussion 


The computations were performed on a cluster machine using 10 to 15 processors. For these transient 
computations at 100,000 ft, because the reaction rates of combustion at such low back pressures are low (compared 
to the rates at sea level), a global time step of 10 ps was used throughout the computations. These global time steps 
used correspond to CFL numbers ranging approximately from 0.1 to unity. The run matrix of the six cases is shown 
in Table 2. The perfectly round case has a unity L/S ratio and zero deformation. The result of the nominal case will 
be discussed first, followed by those of the six out-of-round nozzle extension cases. 


Table 2 Run matrix 


Case 

Description 

L/S ratio 

Deformation, in 

Long axis 

Ref. 

Nominal 

perfectly round 

1.0000 

±0.00 

“ 

19, 30 

1 

slightly out-of-round 

1.0086 

±0.25 

z 

this work 

2 

slightly out-of-round 

1.0086 

±0.25 

y 

this work 

3 

more out-of-round 

1.0346 

±1.00 

Z 

this work 

4 

more out-of-round 

1.0346 

±1.00 

y 

this work 

5 

significantly out-of-round 

1.4400 

±11.6 

Z 

this work 

6 

significantly out-of-round 

1.4400 

±11.6 

y 

this work 


A. The Nominal Case 

The result of the perfectly round case is discussed in here first to establish a baseline that leads to the 
presentation of the out-of-round cases, albeit at a different angle from that of Ref. [19]. That is, instead of showing 
computed total side load, F yZi versus time, computed histories of lateral forces F y and F z are shown in Fig. 4, with the 
timeline of the major side load physics indicated on the right-hand-side of the figure. As indicated by the major 
physics timeline shown in Fig. 4, core jet flow came out as a subsonic core jet after the start command. At 0.4 s into 
the transient, a Mach disk developed and the core jet transformed into a Mach disk flow. Since this initial Mach disk 
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flow is separated from the throat down, the combination of this flow separation pattern and the associated shock 
structure of the Mach disk flow is named the free-shock separation after Nave and Coffey [1], as the supersonic jet 
stemming from the triple point is flowing freely away from the wall. As this Mach disk flow advances downstream, 
the size of the Mach disk grows as the nozzle flowing area increases. At about 0.55 s into the engine start, the initial 


TEG flow begins to 
emerge from the TEG 
exit ring into the nozzle 
extension. Since the 
TEG inlet duct is 
located at zero degree, 
the initial TEG flow is 
asymmetric. The TEG 
flow distribution 

imbalance inevitably 
draws the supersonic jet 
towards the TEM inlet 
duct side of the nozzle 
at around 0.82 s. 



forming an asymmetric 

Fig. 4 Computed side load history during startup for the perfectly round case. 

Mach disk flow with a 


slanted separation line. This phenomenon is unique to this TEM configuration during startup and is named as the 
“asymmetric TEG pumping” in Ref. [19]. As a result, this phenomenon helps the asymmetric FSS Mach disk flow 
transiting into an asymmetric RSS Mach disk flow, generating the first local peak side load at about 0.865s. In a 
RSS flow pattern, first observed and reported during J-2 and J-2S tests, the entire supersonic jet is attached or 
restricted to the nozzle wall [1], Unlike the FSS flow pattern, which has only one separation line, the RSS flow 
pattern has at least one separation line and one reattachment line, and sometimes multiple separation lines and 
reattachment lines. 
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At this time, the rear shock stem is fast approaching the TEG exit ring that composed of a small “base”, or a 
back-step formation between the nozzle contour and the TEG exit ring. As one would expect, when the rear shock 
stem jumps over the base and meets with the flowing TEG flow, a significant disturbance is imposed onto the Mach 
disk flow. By 0.929 s into the startup transient, the rear shock stem, or the separation line, had jumped over the TEG 
exit ring and generated the peak side load of 2114 N. This separation line jump phenomenon was first reported by 
Watanabe [3] as a “separation point jump” during LE-7A engine tests, while Wang [8] captured it later 
computationally in a J-2X transient startup study. Since the separation line jump always produces a peak side load, it 
is considered as the critical side load physics for round, film cooled nozzles. Tomita, et al noted that higher base 
step-height exaggerates the effect of separation line jump [31]. The separation line then moves backup over the TEG 
exit, then coming down again for a second separation line jump at around 0.945 s, with a lower associated side load 
and hence is often ignored. At 1.01 s, the RSS Mach disk flow transitioned back to a FSS Mach disk flow, 
generating another lower local peak side load. This FSS Mach disk flow oscillated several more times in the nozzle 
and eventually left the nozzle at around 1 .49 s. 

While the biggest radius in Fig. 4 indicating when the critical separation line jump occurs. Figure 4 also provides 
a physical aspect that does not show up in the conventional total side load history plot. That is, it shows the lateral 
forces are moving tangentially as the flow moves downstream, or as time accumulates. The tangential motion of the 
lateral forces is only temporarily disrupted by major side load physics such as the separation line jump. After the 
disruption, the tangential motion recovers and resumes its movement along the axial direction. This tangential 
motion of the lateral forces was captured and reported in the SSME startup at sea level study [7], and captured again 
in this study even at 100,000 ft with a shorter flow residence time. The projection of the time-varying locus of the 
lateral forces onto the plane of time zero shown in Fig. 4 gives a conventional, two-dimensional polar plot. 
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Before entering the discussion of the out-of-round cases, a theory that was gradually developed in this study and 
a rationale of using circumferential distribution to demonstrate that theory is best presented first. That is, any 
deformation of the thruster geometry creates a disturbance that makes the flow unstable, or more asymmetrical; On 
the other hand, any deformation also produces a flow directionality that guides the flow in certain favorable 
direction and stabilizes the flow. The balance of these two counter-acting mechanisms decides whether the flow 
becomes more stable, or 


more unstable, which also 
depends on the degree of 
out-of-roundness and the 
design of the nozzle. One 
way to tell which way the 
flow might be going is 
through the interrogation 
of circumferential pressure 
profiles. As shown 
previously in Fig. 2, a 
pressure monitoring ring is 
strategically placed 

between the TEM torus 
and the nozzle extension 

Fig. 5 Computed circumferential pressure profiles at different time slices 
exit plane, such that the during startup for the perfectly round nozzle extension. 

pressure disturbance caused by the geometrical deformation and the pressure stabilization evolved from the flow 
directionality can be sensed. However, it is surmised that when the major side load physics occur, the combined 
pressure disturbances from the side load physics and the geometrical deformation would overshadow that of the 
pressure stabilization owing to the flow directionality. We therefore present the computed circumferential pressure 
distribution at selected time slices that less affected by the major side load events. Three discrete time slices were 
selected for this purpose: 0.65 s, 0.90 s, and 1.80 s, as shown in Fig. 5 for the circumferential pressure profiles of the 
nominal case. By examining the computed results, 0.65 s is a time when the TEG flow fully covers the monitoring 
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ring, yet without any influence from the Mach disk flow, hence it is a good time when the effect of deformation on 
TEG flow can be felt. The time slice of 0.90 s is a snapshot after the shock transition and before the separation line 
jump, hence it is probably a less disturbed time before the disturbance of the critical side load physics arrive. The 
time slice of 1.80 s is about 0.3 s after the nozzle extension flowing full. It is a time that all major side load physics 
have disappeared, and the pressure profile should reflect nothing but the effect of flow directionality, if there is any. 
However, the pressure variations at 0.65 s and 0.90 s in Fig. 5 are just reflections of the effect of lateral forces 
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moving tangentially and showing no sign of flow directionality. The pressure profile at 1.80 s in Fig. 5 is largely 
uniform and showing no signs of flow directionality either. These results are not surprising since there is no flow 
directionality in a round thruster. 


B. Slightly Out-of-round Nozzle Extension Cases 

Figure 6 shows the computed side load histories and the timelines of the major side load physics for the slightly 
out-of-round nozzle extension cases. As shown in Table 2 and Fig. 3, Case 1 is elongated in the z-direction, while 
Case 2 is elongated in the y-direction. As a reminder, these cases represent nozzle extensions permanently deformed 
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due to the internal stresses. The elongation happened in both z- and y-directions since the nozzle extension is thin 
and the material is flexible. Also plotted in Fig. 6 is the main combustion chamber inlet pressure history. It can be 
seen that the pressure ramps up very slowly from about 0.8 s, or spin start, then it ramps up rapidly after about 1,47 
s, or GG start. Since the nozzle extension flows full after about 1.5 s, it can be said that all the major side load 
physics occur during the spin start period. This is because the Mach disk flow is practically evacuated out of the 
thruster by the low ambient pressure at 100,000 ft, with the ambient serves as a vacuum pump. Since the out-of- 
roundness of the nozzle extension is slight, most of the major side load physics that occurred in the nominal case 
also happened in here, although notably missing is the asymmetric TEG pumping. It is speculated that although the 
out-of-roundness is small, the added small directionality is enough to stabilize both the TEG flow and the core flow 
temporarily, thus suppressed the early asymmetric TEG pumping action. Once that short time span is passed, and at 
about 0.83 s into the transient, the Mach disk is big enough such that the supersonic jet feels the slight asymmetry of 
the nozzle extension and becomes more asymmetric (than that of the perfectly round nozzle extension), resulting in 
a more asymmetric FSS-to-Partial RSS (PRSS) transition, and even more unstable separation line jump. The results 
are increased peak side loads of more than 57% and 60%, respectively for Case 1 and Case 2, comparing to that of 
the perfectly round case. Also, due to the slight out-of-roundness of the nozzle extension, the supersonic jet had a 
hard time attach to the wall for the full RSS flow pattern. That resulted in FSS-to-PRSS and PRSS-to-FSS 
transitions. PRSS is a Mach disk flow separation pattern between that of a FSS and a RSS Mach disk flow 
separation patterns. That is, the supersonic jet is only partially attached to the nozzle extension wall. PRSS was first 
captured computationally in the unsteady simulations of LE-7, LE-7A, and CTP50-R5-L nozzles by Wonezawa, et al 
[32] and later captured by Wang [7] in the transient startup computation of an SSME nozzle as well. 
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Figure 7 shows the computed circumferential pressure distributions at selected time slices during startup 
transient for the slightly out-of-round nozzle extension. At 0.65 s, when the pressure profile could show effect of the 
slight flow directionality on the TEG flow, the slight elongation in the z-direction case gives a minimum pressure 
near 90 degree, while the slight elongation in the y-axis gives a minimum pressure near -90 degree. If we reason that 
higher local mass flow results in lower local pressure, then one can see that at 0.65 s, more film coolant flow should 
show up at near 90 degree and -90 degree locations for the slight elongation in the z-direction case, while more film 
coolant flows should prefer 


near the 0 degree and 180 
degree locales for the slight 
elongation in the y-direction 
case. However, the local 
pressure minimums only 
confirmed one quarter of the 
confirmation. Hence, the local 
flow directionality due to the 
slightly out-of-roundness was 
enough to eliminate the TEG 
flow pumping, but not enough 
to overcome the instability 
caused by the slight 
deformation. At 0.90 s, which 
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Fig. 7 Computed circumferential pressure profiles at different time slices 
during startup for the slightly ovalized nozzle extension. 


is just before the separation line jump occurring, the pressure profiles for both cases appear in phase and have local 
pressure minimum at about 40 and -120 degree, and 10 and -90, for the long z-axis and long y-axis cases, 
respectively. Hence, there is no sign of flow directionality at 0.90 s. At 1.80 s, the pressure profiles for both cases 
appear to be near uniform. So, there is no sign of flow directionality at 1.80 s, either. 


C. More Out-of-round Nozzle Extension Cases 
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The more out-of-round nozzle extension cases have a slightly larger L/S ratio of 1.0346 or a ±1 in. physical 
deformation. Such an out-of-roundness is also a likely scenario when there are internal stresses built up through 
many hot-fire tests. Figure 8 shows a comparison of the computed side load histories during startup transient for the 
more out-of-round cases, ovalized and elongated along z- and y-directions, or Case 3 and Case 4. The major side 
load physics depicted in Fig. 8, are similar to those described in Fig. 6. And again, we have PRSS occurring after the 
shock transition due to the out-of-roundness of the nozzle extension. One difference in major side load physics 
worth noting is the long z-axis case has two separation line jumps, whereas the long y-axis case has only one 
separation line jump. It is speculated that since the TEG inlet duct is located at 0 degree, or along the y-axis, and 
intuitively there should be more film coolant flowing under the TEG inlet duct side in the beginning of the TEG 
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Fig. 8 Computed side load history during startup for the more out-of-round Cases 3 and 4. 

emergence, because the TEG inlet duct is closer to the wall for the long axis case than that of the long y-axis case.. 
The long y-axis case is therefore more stable than the long z-axis case, hence it had only one separation line jump. 

Figure 9 shows the comparison of the computed circumferential pressure distributions for Case 3 and Case 4. At 
0.65 s into the transient, this time the peak film coolant flow appears to be near 0 degree location for the long z-axis 
case, while that of the long y-axis case appears to be near the -90 degree location, still no evidence of flow 
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directionality. At 0.91 s into the startup transient, just before the separation line jump, the pressure profiles for the 
two cases appear to be out-of-phase. And, the minimum pressure, or the high flow, occurred at 180 degree location 
for the long z-axis case, while that of the long y-axis case happened near the 0 degree location. Again, there is no 
strong sign of flow directionality for the two cases at 0.91 s. At 1.80 s, however, after the nozzle extension flowing 
full, the effect of higher degree of out-of-roundness is starting to show. That is, there are two local pressure 
minimums occurring for the long z-axis case, at near 90 degrees and -90 degrees, respectively. In the mean time, 
there are two local pressure minimums occurring for the long y-axis case as well, at near 0 degrees and 180 degrees. 
The flow directionality definitely shows for the 1.80 s pressure profiles, since the effect of transient flow phenomena 
are long gone. I can also be said 


that the effect of flow 
directionality of the more out-of- 
round cases is slightly bigger 
than that of the slightly out of- 
round cases during the 
separation line jump, just it was 
overshadowed slightly by the 
effect of deformation. Hence, 
the peak slide loads for Case 3 
and Case 4 only increased by 
50% and 55%, respectively; 
slightly less than those of 57% 



and 60% for Case 1 and Case2, 
respectively. 


Fig. 9 Computed circumferential pressure profiles at different time 
slices during startup for the more ovalized nozzle extension. 


D. Significantly Out-of-round Nozzle Extension Cases 

The results of the last two cases, or Case 5 and Case 6, are the hypothetical, significantly out-of-round nozzle 
extension cases. With a long-axis-to-short-axis ratio of 1.4400, and a maximum physical deformation of ±1 1.6 in at 
the nozzle extension exit plane, it is anticipated that the flow stability improves with the increased flow 
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directionality. As a result, the computed peak side loads should be lower than those of the four afore -mentioned 
cases. Figure 10 shows the computed side load histories during startup transient for the significantly out-of-round 
cases. Indeed, the computed peak side loads are now +28% and -18% from those of the nominal cases, and much 
lower than those of the more out-of-round cases. In addition, due to the significant ovalization of the nozzle 
extensions, instead of the having the major side load physics consisting FSS-to-PRSS transition, separation line 
jump, and PRSS-to-FSS transition; we now have major side load events comprising FSS-to-PRSS transition, 
separation line jump, PRSS-to-FSS transition, FSS-to-PRSS transition, and PRSS-to-FSS transition. That is, two 
more shock transitions before the nozzle extension is flowing full. Furthermore, in the long z-axis case, we have 
three separation line jumps, while there is only one separation line jump for the long y-axis case. Again, this is 
because the TEG inlet duct is closer to wall in the long z-axis case, while the TEG inlet duct is farther away from the 
wall. As a result. Case 6 is more stable than Case 5. For the same reason, the peak side load of Case 6 is 
significantly less than that of Case 5. 

The reason of the peak side load of Case 6 is lower than that of Case 5 can also be explained by the computed 
circumferential pressure distributions at selected time slices of 0.65 s, 0.92 s, and 1.80 s, as show in Fig. 1 1. At 0.65 
s, when there is no interaction between the TEG flow and the Mach disk flow, the pressure profiles are solely 
coming from the contribution of the TEG flow. It can be seen that there are two local pressure minimums located at 
90 degree and -90 degree for the long z-axis case, and two local pressure minimums located at 0 degree and 180 
degree for the long-y-axis case. This exhibition of clear flow directionality can be visualized from the cross- 
sectional view of Case 5 and Case 6 in Fig. 3, where the least flow resistances to the film coolant for the long z-axis 
case are the 90 degree and -90 degree locations. Similarly, the least flow resistances to the TEG for the long y-axis 
case are the 0 degree and 180 degree locations. Recall that similar phenomenon was not revealed until at 1.80 s for 
the more out-of-round cases in Fig. 9. In the significantly out-of-round cases, this phenomenon, or strong flow 
directionality due to significant ovalization of the nozzle extension, just happened much earlier, at 0.65 s. The effect 
of flow directionality occurs earlier as the degree of out-of-roundness increases. 


21 



At 0.92 s into the startup transient, just before the occurrence of the peak side load physics - the jump of the 
separation line, one can see that the long y-axis case has two local pressure minimums, located at approximately 0 
degree and 180 degree, showing strong influence from the flow directionality; while the long z-axis case has only 
one local pressure minimum, located at around 120 to 180 degrees, showing no sign of flow directionality. That 
means that at this time, the flow of the long y-axis case is more symmetric and more stable, than that of the long z- 
axis case, assisted by more flow directionality, just before the separation line jump. As a result, the peak side load 
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Fig. 10 Computed side load history during startup for the significantly out-of-round cases (L/S = 1.4400). 


due to separation line jump of the long y-axis is much lower than that of the long z-axis case; in fact, 18% lower 
than that of the nominal case. Finally, at 1.80 s into the startup transient, when the nozzle extension is flowing full, 
the effect of flow directionality manifested itself strongly. The flow at the long z-axis case clearly favors the 90 
degree and -90 degree directions, while that of the long y-axis case clearly prefers the 0 degree and 180 degree 
directions. The flows at the two cases are now symmetric. 
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The major flow physics of case 6 are now presented as Mach number contours at selected time slices, as shown 


in Fig. 12. For each time slice, two pictures are shown. The picture on the left-hand-side is the Mach number 
contours on the xy-plane. 


while that on the right- 
hand-side is the Mach 
number contours on the xz- 
plane. The xy-plane views 
the long-axis side, or wider 
dimension of the ovalized 
nozzle, and the xz-plane 
views the short-axis side, or 
narrower dimension of the 
ovalized nozzle. During the 
initial startup transient, the 
FSS Mach disk flow starts 
at about 0.39 s, while the 



Fig. 11 Computed circumferential pressure profiles at different time slices 
during startup for the significantly out-of-round nozzle extension. 


TEG emerges at about 0.55 

s. At 0.61 s into the startup transient, the Mach disk of the FSS Mach disk flow is not large enough to have the 
supersonic jet affecting the TEG flow, and vice versa. Until about 0.80 s, the Mach disk is large enough and the 
TEG flow and the supersonic jet begins to attracting each other. At 0.848 s, the supersonic jet and the TEG flow 
merge and PRSS Mach disk flow is formed. The Mach number contours at 0.865 s in Fig. 12 shows a typical PRSS 
Mach disk flow where the supersonic jet is partially attached to the nozzle extension and the reattachment is not 
continuous. At 0.925 s, about 0.03 s before the occurrence of the separation line jump, the top separation line is 
slightly wavy, and the bottom reattachment line is continuous, hence this is a rare RSS Mach disk flow for this case. 
The separation line jump occurred at 0.928 s and it triggered a peak side load of 1738 N. At 0.93 s, or 0.02 s after 
the critical side load, the Mach disk flow is quite disturbed, and there are three discontinuous separation lines and 
two discontinuous reattachment lines. Note the shock stem (foot) jumped over the TEG exit ring on the xy-plane, 
and lagged behind the exit ring on the xz-plane. As a result, there is a partial separation line hidden behind the TEM 
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torus on the xy-plane. That partial separation line is the residual of the top, wavy separation line at 0.925 s that had 
mostly jumped over the base formation of the TEG inlet and disappeared thereafter - or the remaining separation 
line of the splitting of the separation line mechanism [19]. Afterwards, sickle-shaped separation line emerged due to 
the deformation of the nozzle extension. That sickle-shaped separation line flow pattern was temporarily disrupted 
to become circular separation line flow pattern after the first PRSS-to-FSS transition (not shown), but resumed the 
formation at 1.095 s and maintained that formation until the nozzle extension was flowing full. The Mach number 
contours at 1.125 s time slice in Fig. 12 shows an initial sickle is growing at the two ends of the short axis. At 1.4 s 
into the startup transient, the sickle grows to its largest and the end of the sickle reaches the end of the nozzle 


extension. 


A comparison of the computed 
peak side loads is summarized in 
Table 3, along with those computed 
in a previous study [19]. In this 
study, comparing to the peak side 
load of the nominal case, those have 
increased 57% and 60% for the 
slightly out-of-round cases, 50% and 
55% for the more out-of-round 
cases, and 28% and -18% for the 
significantly out-of-round cases. In 
previous study, comparing to the 
peak side load of the perfectly round 
case, that has increased 26% for the 
slightly out-of-round case, 55% for 
the more out-of-round case, and only 
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round case. In general, the overall 


peak side loads culminated at the slightly out-of-round case in this study, and climaxed at the more out-of-round 
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case in previous study. This is caused by the difference in the way the thrusters were ovalized, such that the flow 
instability overpowers the flow directionality at the slightly out-of-round case in this study, while that happens at the 
more out-of-round case in previous study. 


Table 3. A comparison of the computed peak side loads, kN 


Case 

Description 

Long axis 

This study 

Previous study 

nominal 

perfectly round 

- 

2114 [19] 

2114 [19] 

1 

slightly out-of-round 

z 

3309 (+57%) 

2668 (+26%) [19] 

2 

slightly out-of-round 

y 

3376 (+60%) 

- 

3 

more out-of-round 

Z 

3175 (+50%) 

3275 (+55%) [19] 

4 

more out-of-round 

y 

3268 (+55%) 

- 

5 

significantly out-of-round 

Z 

2715 (+28%) 

2171 (+2.7%) [19] 

6 

significantly out-of-round 

y 

1738 (-18%) 

- 


Finally, the discussion would be incomplete without mentioning the study performed by Ostlund and Bigert [9]. 
They studied several non-round, polygon nozzles. Those polygon nozzles were not deformed nozzles, but purposely 
designed as faceted nozzles in order to reduce nozzle side loads. One of the side load reduction mechanisms 
pertinent to the current study is that the polygon corners acted as a kind of structure-breaker, splitting the flow 
separation pattern in the azimuthal direction. That mechanism, splitting the flow separation pattern, is essentially 
equivalent to the splitting of the separation line mechanism reported in previous study [19] for the significantly out- 
of-round case, which also resulted in side load reduction. That agreement in side load reduction mechanism not only 
provided experimental evidence that the computed flow pattern [ 19] was credible, but also the mechanism itself is in 
agreement with the general theory developed in this study for out-of-round nozzles. In other words, splitting of the 
separation line [19] is the result of the increased flow directionality in the elongated side, made possible by the 
significant our-of-roundness over the entire thrust chamber [19] or over the nozzle extension as discussed in this 
study. Indeed, the faceted polygon nozzles provided increased flow directionality in the polygon corners, thereby 
stabilizing the inherently unstable flow in the original perfectly round nozzle, achieving nozzle side load reduction. 
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In summary, the flow in the perfectly round case is inherently unstable, since there is no geometrical deformation 
to provide any flow directionality to stabilize the flow. In the mean time, any out-of-roundness, whether it is 
deformed in the entire thrust chamber, or the nozzle extension only, boosts the flow destabilization since any 
geometrical out-of-roundness exacerbates the flow asymmetry. On the other hand, any out-of-roundness also 
provides a flow directionality, which eventually stabilizes the flow. That is to say, the deformation generates two 
counteracting mechanisms, one destabilizes the flow that causes the side load goes higher, another stabilizes the 
flow that reduces the side load. In this effort, we see the flow destabilizing effect is stronger in the slightly out-of- 
round cases and the computed circumferential pressure profiles do not show discernible effect of flow directionality; 
In the more out-of-round cases, the flow stabilizing effect gains strength and we see obvious flow directionality in 
the computed circumferential pressure distributions at time slice of 1.80 s; In the significantly out-of-round cases, 
the flow stabilizing effect grows strength and signs of flow directionality are noticeable in the computed 
circumferential pressure profiles in all three time slices. 

VI. Conclusion 

Three-dimensional numerical computations were performed to determine the effect of out-of-round nozzle 
extension on nozzle side loads during engine start transient for a film cooled nozzle. The out-of-roundness was 
achieved with ovalization of the nozzle extension. Six cases were performed with three different degrees of 
ovalization on y and z axes: one slightly ovalized, one more ovalized and one significantly ovalized. It was found 
that the separation line jump is the critical side load physics for all cases. Comparing to the peak side load of the 
perfectly round nozzle, those of the slightly out-of-round cases increased by 57% and 60%, those of the more out-of- 
round case increased by 50% and 55%, and those of the significantly out-of-round cases increased by 28% and - 
18%, for the elongation on the z- and y-axis, respectively. For the significantly ovalized nozzle extension, the 
difference between the peak side loads of the two elongation cases comes from the location of the single TEG inlet 
duct, which resides on the y-axis. 

Based on the observations from the results of this study and two other references, a general theory on the impact 
of geometrical deformation on the transient nozzle side loads is developed. That is, the deformation increases the 
flow instability in the nozzle, which increases the flow asymmetry and the nozzle side load. On the other hand, the 
deformation also increases the flow directionality, which decreases the flow instability and the nozzle side load. The 
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final outcome between the two counteracting mechanisms depends not only on the degree of ovalization, but also on 
the specific type of deformations and design of nozzle. Under the circumstances of this effort, the flow instability 
overtakes the flow directionality in the slightly ovalized nozzle extension cases, while sign of flow directionality 
developed after nozzle extension flowing full in the more ovalized nozzle extension cases, and the flow 
directionality overcomes the flow instability in the significantly out-of-round nozzle extension cases. 
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What is nozzle side load? 


• Nozzle side load is a lateral force occurring naturally 
during transient startup and shutdown of rocket engines 

• Nozzle side loads are detrimental to the engine 
components such as actuator, cooling tubes, propellant 
feed lines, etc. for almost all rocket engines during 
development 

• For regeneratively cooled engines, such as SSME, the 
origins of the side loads have been identified as 
combustion wave, shock transitions, and shock breathing 
at the nozzle exit 

• For film cooled engines, such as LE-7A and J-2X, the 
origins of side loads have been identified as shock 
transitions, and flow separation line jump over the turbine 
exhaust gas (TEG) flow 
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J-2X Engine 


• J-2X engine, the Ares I upper 
stage engine under 
development, is an evolved 
variation of two historic 
predecessors: the J-2 and J-2S 
engines. 

• It has a regeneratively cooled 
combustion chamber and 
nozzle, along with a film-cooled 
nozzle extension. 
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Nozzle Out-of-Roundness 


• Large, relatively light weight, liquid rocket engine 
nozzles, are probably never truly round. 

• The causes of nozzle out-of-roundness could be, but are 
not limited to, long term internal stresses, especially high 
thermal load during hot-fire test, such as creep. 


Objective 


• To computationally investigate the effect of out-of-round 
nozzle extension on the nozzle side load during the J-2X 
engine startup transient. 
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Computational Methodology 


• Multidisciplinary computational methodology: 


■ UNIC is a unstructured-grid, pressure-based, reacting 
flow, CFD & heat transfer solver 

■ Engine system modeling for transient inlet properties (to 
simulate hot-firing tests) 

■ Thermal modeling of wall temperatures for combustion 
chamber, nozzle, and nozzle extension (to simulate hot- 
firing tests) 
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Benchmark with Regeneratively Cooled SSME Nozzle 
during Sea Level Startup 
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Benchmark with Regeneratively Cooled SSME Nozzle 


During Sea Level Startup 
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Com parison of J-2X Sea Level Results w ith those of LE-7A 

• Watanabe’s LE-7A hot-fire Tests at sea 
level (Watanabe, Yasuhide, et al. , “LE- 
7A Engine Nozzle Problems during the 
Transient Operations,” AIAA Paper 
2002-3841). 

■ Measured the transient side loads of LE-7A 
engine with and without the nozzle extension 

■ Reported that peak side load was caused by 
the separation line jump 

■ Found out the side load of the stub nozzle 
was reduced drastically. 




Comparison of Computed J-2X (Nozzlette) with 
those of LE-7A Hot-Fire Test 


Comparison of the Sea Level Peak Side Loads 
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A Grid Layout of the J-2X Nozzle 


Engine 

Grid size 

Circumferential 
division number 

SSME 

1,275,120 

72 

J-2X 

Nozzlette 

2.058192 

72 

J-2XTEM 

4,421,166 

120 

J-2XTEM 
& Duct 

20,222,021 

240 
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TEM torus 


Turbine Exhaust 
Manifold (TEM) 
inlet duct 


MCC 


Nozzle 


Nozzle 

extension 


Pressure- 

monitoring 

ring 


A separate steady-state mainstage analysis 
with this TEM configuration indicated it 
produces a 5% higher, than uniform, mass 
flow into the nozzle extension in the region 
under the TEM inlet duct. 






Cross-Sectional View of the Six Ovalized Thrusters 
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Run Matrix 


& 


Case 

Description 

L/S ratio 

Deformation, in 

Long axis 

Ref. 

baseline 

perfectly round 

1.0000 

±0.00 

- 

19,30 

1 slightly out-of-round 1.0086 ±0.25 Z this work 

2 

slightly out-of-round 

1.0086 

±0.25 

Y 

this work 

3 more out-of-round 

1.0346 ±1.00 Z 

this work 

4 

more out-of-round 

1.0346 

±1.00 

Y 

this work 

5 significantly out-of- 1.4400 ±11.6 Z this work 

round 

6 

significantly out-of- 
round 

1.4400 

±11.6 

Y 

this work 


14 











Transient Startup Inlet Flow Properties of MCC and TEG 
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Computed Side Load History for Nominal Case 
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Computed Circumferential pressure Histories 
for Slightly Overali/ed Nozzle Extensions 
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Computed Side Load Histories for 
More Ovalized Nozzle Extensions 
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Computed circumferential pressure Histories 
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Computed Side Load Histories for 
Significantly Ovalized Nozzle Extension 
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Computed circumferential pressure Histories 
for Significantly Overalized Nozzle Extensions 
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Computed Mach Number Contours 
for Significantly Overalized Nozzle Extensions 
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A Comparison of the Computed Local Peak Side Loads 



Case 

Peak side loads, kN 


Description 

Long axis 

This study 

Previous study 

Nominal 

perfectly round 

- 

2114 [19] 

2114 [19] 

1 

slightly ovallized 

z 

3309 (+57%) 

2668 (+26%) [19] 

2 

slightly ovalized 

y 

3376 (+60%) 

- 

3 

more ovalized 

z 

3175 (+50%) 

3275 (+55%) [19] 

4 

more ovalized 

y 

3268 (+55%) 

- 

5 

significantly 

ovalized 

z 

2715 (+28%) 

2171 (+2.7%) [19] 

6 

significantly 

y 

1738 (-18%) 

- 


ovalized 
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Conclusions 


• Six out-of-round nozzle extension cases were computed with three different 
degrees of ovalization on y and z axes: a slightly ovalized, a more ovalized, 
and a significantly ovalized. 

■ Separation line jump is the critical side load physics 

■ Comparing to the peak side load of the perfectly round thruster, those of the 
slightly ovalized cases increased by 57% and 60%, those of more ovalized 
cases increased by 50% and 55%, and those of the significanly ovalized 
cases increased by 28% and -18% for the elongation along the z and y 
axes, respectively. 

■ For the significantly ovalized cases, the difference between the peak side 
loads comes from the location of the single TEG inlet duct. 

• A general theory on the impact of geometry deformation on the transient peak 
side loads is developed: 

• Under the assumptions of this effort, the flow instability overtakes the flow 
directionality in the slight ovalized cases, while sign of flow directionality 
developed after thruster flowing full in the more ovalized cases, and the 
flow directionality overcomes flow instability in the significantly ovalized 
cases. 
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